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We revisit the origin of structures problem of standard Friedmann-Robertson- 
Walker cosmology to point out an unjustified approximation in the prevalent analysis. 
We follow common procedures in statistical mechanics to revise the issue without 
the disputed approximation. Our conclusions contradict the current wisdom and 
reveal and unexpected scenario for the origin of primordial cosmological structures. 
We show that standard physics operating in the cosmic plasma during the radiation 
dominated expansion of the universe produce at the time of decoupling scale invariant 
density anisotropics over cosmologically large comoving volumes. Scale invar iance is 
shown to be a direct consequence of the causality constrains imposed by the short 
FRW comoving horizon at decoupling, which strongly suppress the power spectrum 
of density fluctuations with cosmologically large comoving wavelength. The global 
amplitude of these cosmological density anisotropies is fixed by the power spectrum in 
comoving modes whose wavelength is shorter than the causal horizon at the time and 
can be comparable to the amplitude of the primordial cosmological inhomogeneities 
imprinted in the cosmic microwave background radiation. 



I. INTRODUCTION. 



The cosmic microwave background radiation (CMBR) is a fossil evidence of the early his- 
tory of the expanding and cooling Friedmann-Robertson- Walker (FRW) universe. Indeed, 
this radiation is a high-resolution picture of the cosmological density field at the instant of 
recombination, when the temperature of the universe was T eq ~ 1 eV and this radiation last 
scattered and decoupled from the cosmic plasma. At that instant the size of the universe 
was 10~ 3 of its current size. The red shifted radiation observed at present is extremely 
homogeneous and isotropic over the sky, with a perfect blackbody spectrum of temperature 



T = 2.725 K ~ 10" 3 eV Q up to tiny fluctuations AT/T < 10" 5 Q, ||. Encoded in the 
statistical correlations of these tiny temperature anisotropies are the statistical features of 
a primordial seed of very small cosmological density anisotropies in the, otherwise, homoge- 
neous and isotropic universe at the instant of decoupling. Current cosmological models work 
within the paradigm that large scale mass structures in the distribution of matter in the 
present universe are the result of gravitational evolution operating since the time of matter 
radiation equality on that primordial seed (for a recent claim of observational evidence of 
this correlation see [^]). 

The most characteristic statistical feature of the primordial seed of density anisotropies 
at decoupling is its scale invariance over cosmologically large comoving scales the vari- 
ance (AMy) 2 of density anisotropies over volumes V of cosmological comoving size L is 
proportional to the area S ~ L 2 of the surface that bounds the considered region, 

(AMy) 2 ~ k 2 p % H-'S. (1) 

The dimensionless numerical factor k ~ 10 -5 measures the global amplitude of the primor- 
dial density anisotropies in terms of the cosmological density at decoupling p eq = ^gf- and 
the comoving causal horizon at that time. Comoving scales of cosmological size are 
much larger than this horizon, 

H- 1 « L < Ho 1 . (2) 

The upper bound Hq 1 ~ lO 4 ^" 1 is the comoving size of the present causal horizon. For 
the sake of simplicity we have chosen the scale factor of the expansion of the universe at 
decoupling a(t eg ) = 1 as reference. 

Scale invariant macroscopic density anisotropies ((TJ) produce a gravitational potential 
independent of the scale L, 

(A0U) 2 = G 2 ^l . n\ (3) 

This feature is observable, through the Wolfe-Sachs effect 0, as the renowned Harrison- 
Zeldovich plateau in the spectrum of temperature anisotropies in the CMBR at angular 
scales larger than the angle subtended on the sphere of last scattering by the comoving 
horizon at decoupling, 
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-A(f) ~k, at 9 > 1°. 
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The amplitude of temperature fluctuations over large angular scales is roughly constant 



(4) 



A. The prevalent analysis. 

The pattern (1 1121) of scale invariant primordial density anisotropies at decoupling over 
comoving volumes of size much larger than the causal horizon at that time has puzzled 
cosmologists since long ago. Theoretical analysis reproduced in textbooks as well as research 
papers |& 0, IllJ seem to conclusively prove that the pattern cannot be reproduced in the 
(otherwise very successful) standard FRW cosmology without violating causal correlations 
over scales beyond the horizon. The argument goes as follows (see, for example, Eq. 9.1, 
9.2 and 9.3 in Q). 

The variance of macroscopic density anisotropies at decoupling over volumes V of co- 
moving size L, which is related to their power spectrum V{k) by the equality (see section 
II.B) 



(AM,) 2 = p% 
is approximated by the expression: 
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Approximation © relies on the suppression of the geometric factor (JZJ) for modes k ^> 
1/L to estimate that the contribution of these modes to macroscopic inhomogeneities (JSJ 
over comoving volumes of size L is negligible and conclude that the dominant contribution 
comes necessarily from density fluctuations in modes k < 1/L with comparable comoving 
wavelength. 
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According to expression (JHJ) the pattern (11121) of scale invariant primordial density 
anisotropies over cosmologically large comoving volumes corresponds univocally to a power 
spectrum at decoupling 

V*(k) ~ (k 2 H~ a ) k, for k < H eq , (8) 

which decreases linearly over the range of modes with cosmologically short comoving mo- 
menta. To check this inference plug (JHJ) into (JHJ) to reproduce the scale invariant pattern 
(HE): (AM v )l ~ pl q V*(k ~ 1/L) V ^ k 2 pl q H-^ V ~ k 2 p 2 eq H~ q S. 

Inference (JHJ) is bewildering within the theoretical framework of standard FRW cosmology 
The power spectrum of density anisotropies that standard physics operating during the 
radiation dominated expansion of the universe can produce at decoupling is known to be 
suppressed much strongly than linearly (JHJ) over the range k <C H eq . The suppression is 
known to be at least quartic [ic| . 

V{k) < V{k ~ H eq ) (J^J , far k<^H eq . (9) 

This feature is forced by the constrain of global mass (energy) conservation and by the causal 
structure of the FRW universe, in which sites separated by cosmologically large comoving 
distances have not ever been in causal contact before decoupling. Hence, (JSJ) applies to the 
power spectrum of density anisotropies generated by any mechanism operating in the cosmic 
baryon-photon plasma before decoupling as long as it does not violate causality. 

A direct evaluation, according to approximation (JSJ, of the variance of cosmological 
density anisotropies at decoupling that correspond to the power spectrum (JOJ) leads to 

(AM,) 2 ~ (AM V J 2 ^. (10) 

According to this estimation, the variance of density anisotropies that standard physics 
produce at decoupling over cosmologically large comoving volumes V of size L ^> H~ q l is 
predicted to be suppressed by a factor H~ q / L with respect to the variance of the anisotropies 
(AM Veq ) 2 that the same physics generate at that time over comoving volumes V eq of the size 
of the causal horizon. The manifest disagreement between the theoretical prediction (JTUJ) and 
the observed scale invariant pattern (0) is known in the literature as the origin of structures 
problem of standard FRW cosmology. 
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In order to solve this problem cosmologists have worked out during the last twenty years 
a quite unexpected cosmological paradigm: a period of exponential inflation when the tem- 
perature of the universe was close to the Planck scale would have supposedly preceded the 
standard radiation dominated expansion of the universe It is claimed that a very early 
period of cosmological inflation would naturally explain the exceptional conditions of the 
universe at the beginning of the standard FRW expansion: namely, its huge size and its 
flatness, as well as its homogeneity and isotropy over cosmologically large comoving scales. 
In addition, the density anisotropics present in the very first instants of the history of the 
universe would be stretched by the inflation to comoving cosmological scales. At decoupling 
these anisotropies would have a characteristic linear power spectrum 

V inflat {k) ~ A k, for k^H eq (11) 

over the range of cosmologically short comoving momenta. According to inference (jHJ) these 
anisotropies would be scale invariant, (AM v ) 2 n f lat ~ p 2 eq Vi n fi at {k ~ 1/L) V ~ p 2 eq A S. The 
global amplitude of the anisotropies left by inflation is fixed by the parameter A. In most 
models of inflation, A is an adjustable parameter that depends on the scale at which inflation 
takes place. This parameter is tuned to the value A ~ k 2 H~ q by additional constrains on 
the inflationary setup. 

The widely accepted conclusion of the analysis that we have briefly reviewed is that 
primordial scale invariant cosmological anisotropies (1 1121) are a pristine signature of inflation, 
which cannot have been altered by anisotropies (jl())l generated by standard physics during 
the radiation dominated expansion of the universe. Inflation is currently regarded by most 
cosmologists as the natural solution to the origin of structures problem of FRW cosmology, 
even though multiple technical problems (mostly related with the need to causally connect 
two temperature scales T eq ~ 1 eV and T P i anck ~ 10 18 GeV separated by 28 orders of 
magnitude) have prevented the inflationary paradigm to get successfully implemented in 
models of particle cosmology. 

B. An unjustified approximation. 

The motivation of this paper is to contend that the above analysis of the origin of struc- 
tures problem of standard FRW cosmology is flawed by an unjustified approximation. We 
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point out that the estimation (JTUj) of the variance of cosmological density anisotropies pro- 
duced at decoupling by standard physics is deceptive and leads to erroneous conclusions. 
The estimation (|TU|) is based on approximation which relies on the suppression of the 
geometric factor (JZj) in modes k 3> 1/L in order to neglect the contribution of these modes 
to the integral expression (JSJ) for the variance of density anisotropies over comoving volumes 
of size L. According to estimation (fTU|) . the suppression of the geometric factor in modes 
k ^> 1/L implies that the largest contribution to the integral expression (jSJ) must come from 
modes k < 1/L. We notice that this argument is invalid when the considered comoving 
volume of integration is cosmo logically large L ^> H~g, because the power spectrum (jHJ) 
of density anisotropies left at decoupling by standard physics is greatly enhanced in modes 
k ~ H eq 3> 1/L and suppressed in modes k < 1/L, 

V(k ~ 1/L) < V{k~ H eq ) f-gLj , for L^H~\ (12) 

In order to get a reliable estimation of the contribution of modes k > H eq to the variance (jSJ) 
of density anisotropies produced at decoupling by standard physics (0) over cosmologically 
large comoving volumes of size L ^> if" 1 , we must compare the suppression of the geometric 
factor ((7j) in these modes against their enhanced power spectrum ()12|). 

Indeed we have proven in [l^ j. Section III, that approximation © does not give a cor- 
rect estimation of the variance (0) of cosmological density anisotropies whenever the power 
spectrum V(k) ~ k n decreases faster than linearly, n > 1, over the range k H eq of 
comoving modes with cosmologically large comoving wavelength. For these power spectra 
the largest contribution to the variance of cosmological density inhomogeneities comes from 
modes k > H eq with comoving wavelength within the horizon, while the contribution from 
modes k < 1/L with cosmologically large comoving wavelength is negligible. Approxima- 
tion (jUJ) carelessly neglects the contribution of modes with comoving wavelength within the 
horizon to the variance of density anisotropies over cosmologically large comoving volumes. 
This forgotten contribution accounts for random density fluctuations at the boundary sur- 
face of the cosmological volume of integration. Once this contribution is taken into account 
the variance (0) of cosmological density anisotropies associated to the power spectra Q is 
shown to be scale invariant, 
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This theoretical prediction mends the flawed estimation ()10[). In next Section II, subsections 
II. B and II. C, we prove and discuss this prediction in clear physical terms. 

In summary, we claim that standard physics operating during the radiation dominated 
expansion of the universe generate at decoupling scale invariant density anisotropies over cos- 
mologically large comoving volumes. The cosmological anisotropies are produced by random 
density fluctuations with comoving wavelength shorter than the comoving causal horizon at 
the time. According to prediction (jlHj) . the pattern (1 1121) of scale invariant primordial cos- 
mological anisotropies is not bewildering in the framework of standard FRW cosmology. 
On the contrary, this pattern is the clear signature of the strong suppression of the power 
spectrum of standard physics in modes with cosmologically large comoving wavelength. 
Scale invariance is dictated by the causality constrains of standard FRW cosmology. This 
conclusion obviously contradicts the current wisdom in cosmology, expressed in the flawed 



Equation ()13j) also proves that the inference (JHJ) of a linear power spectrum of primordial 
cosmological density anisotropies at the time of decoupling is absolutely unjustified. Scale 
invariance (1 1121) of primordial cosmological density anisotropies at decopling only implies 
that their power spectrum decreases faster than linearly over the range of cosmologically 
short comoving momenta. The power spectrum @ fulfills this constrain. 



The global amplitude of the scale invariant cosmological density anisotropies (j!3)l pro- 
duced at decoupling by standard physics operating during the radiation dominated expan- 
sion of the universe is fixed by their power spectrum V(k) in modes k > H eq with comoving 
wavelength shorter than the causal horizon at that time. Hence, their amplitude is propor- 
tional to the amplitude (AMy e? ) 2 ~ p 2 eq V{k ~ H eq ) V eq of the inhomogeneities that the same 
physics produce over comoving volumes V eq of the size of this causal horizon. A comparison 
of ()13|) with (^) implies that it is necessary 
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C. A revised analysis. 
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for standard physics operating during the radiation dominated expansion of the universe to 
produce at decoupling scale invariant cosmological density anisotropies with global ampli- 
tude comparable to that of the primordial scale invariant density inhomogeneities imprinted 
in the CMBR. This conclusion reveals an unexplored, but feasible, scenario for the origin 
of primordial cosmological structures within the framework of FRW cosmology Q|. This 
scenario had been previously disregarded in the literature, in favour of the inflationary 
paradigm, due to the flawed estimation (JTUJ). 

During the radiation dominated expansion of the universe the cosmic plasma is a dense 
and opaque soup of coupled baryons and photons in thermal equilibrium. Density (pressure) 
waves can propagate in this fluid with sound velocity c s ~ (in natural units). These waves 
homogenize the cosmic density field, eroding all random density fluctuations with comoving 
wavelength much shorter than the sound causal horizon at the time. Obviously, causality 
constrains prevents this mechanism to operate over comoving distances larger than the sound 
causal horizon. Hence, the power spectrum of the density anisotropies left by standard 
physics in the cosmic plasma at decoupling is peaked at k ~ H eq and suppressed V(k) < 
V{k ~ H eq ) f^ 21 ) over modes k 3> H eq by the action of these non-linear wave phenomena. 
We have discussed above Q that the power spectrum of the density anisotropies left by 
standard physics is also suppressed over modes k -C H eq because of causal constrains. Such 
power spectrum corresponds to a pattern of macroscopic density anisotropies (AMy) 2 ~ 
(AMv e J 2 (y~) over volumes of comoving size L < H~ q , which turns scale invariant ()13|) 
over volumes of comoving size L 3> H~ q . We will show in next Section II, subsection II. E, 
that the amplitude (fTljl is roughly given by 

(AM Ve f ~ 2([ P (f)] 2 - pi) vl = 2 P % ( [ ^v(k) | \ rns 



where p(x) is the cosmological density field at decoupling. The dimensionless factor can be 

d 3 k 
(2tt) 3 



estimated by / 4^V{k) ~ V(k ~ H eq ) H* q . Thus, Eq. (JUJ) is satisfied if 



V(k ~ H eq ) H% — — 10" 10 . (16) 

The dimensionless factor / j^is'Pik) that fixes (|15j) the global amplitude of the cos- 
mological density anisotropies at decoupling (fE?j) ia a measure also of the entropy of the 
fluctuations. If the radiation dominated expansion of the universe is adiabatic and the en- 
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tropy of the density fluctuations is conserved, then the global amplitude of the cosmological 
density anisotropies at decoupling is set by the dynamics of the universe at a much earlier 
stage. We discuss condition (|THj) in further detail below in Section II, subsection II. E. 

The implications of this revised analysis are far reaching. The standard FRW cosmolog- 
ical framework that we are discussing here assumes, implicitly, that the exceptional initial 
conditions at the beginning of the radiation dominated expansion of the universe (namely, its 
huge size and flatness, as well as its homogeneity and isotropy over cosmologically large co- 
moving scales) have been set by some unspecified mechanism (e.g. inflation) at some earlier 
stage. According to the prevalent analysis of the origin of primordial cosmolgical struc- 
tures, see subsection LA, the mechanism that fixes the initial homogeneous and isotropic 
conditions of the universe needs to be tuned to seed at once density fluctuations with cosmo- 
logically large comoving wavelength and the right amplitude, see Eq. (jl 1)1 . The difficulties 
to get the mechanism appropriately tuned have an obvious origin: cosmological homogene- 
ity and isotropy means, by definition, that all density fluctuations with cosmologically large 
comoving wavelength should be negligibly small. We realize now that, according to Eq. 

(USED, 

• The primordial cosmological density anisotropies (1 11211 have not to be necessarily 
seeded by the same mechanism that fixes the initial homogeneous and isotropic 
conditions at the beginning of the standard FRW expansion of the universe. The 
anisotropies may be seeded later on, during the radiation dominated expansion of 
the universe. Thus, we should consider as feasible models of inflation (or any other 
mechanism) that can fix the initial size, flatness, homogeneity and isotropy of the 
FRW universe, and leaves negligibly small density fluctuations with cosmologically 
large wavelength. Such models should be easier to implement in models of particle 
cosmology. The price to pay is obvious: the physics of such a mechanism would not 
be readily testable on the spectrum of temperature anisotropies in the CMBR. 

• Scale invariant primordial density anisotropies (jll2j) over cosmologically large comov- 
ing volumes are not necessarily a pristine signature of relic inhomogeneities left over 
by the mechanism that fixed the initial conditions at the beginning of the radiation 
dominated expansion of the universe (e.g. inflation). Even in the case that relic scale 
invariant cosmological inhomogeneities were left by that mechanism at the beginning 



10 



of the radiation dominated expansion, standard physics operating later on produce 
scale invariant cosmological anisotropics that are overlaid upon them by the time of 
decoupling. 



II. SCALE INVARIANT DENSITY ANISOTROPIES. 

The main aim of this paper is to prove that the characteristic scale invariance of the 
primordial cosmological density anisotropies at decoupling ()l|2j) is an unequivocal signature 
of the standard FRW cosmology. In the standard cosmology the causal comoving horizon 
grows monotonically with cosmic time. Comoving scales of cosmological size are much larger 
than the comoving causal horizon at decoupling. The short comoving causal horizon pre- 
vents the generation of density fluctuations with cosmologically large comoving wavelength 
during the radiation dominated expansion of the universe (jUJ). As a consequence the domi- 
nant contribution to density anisotropies at decoupling over cosmologically large comoving 
volumes comes from random density fluctuations with comoving wavelength shorter than 
the horizon. This contribution corresponds to random density fluctuations at the boundary 
surface of the volume of integration and, hence, the variance of these anisotropies grows 
with the area of the boundary surface. Below we compare the geometric dependence of the 
variance of the primordial cosmological density anisotropies on the boundary area with 
the dependence of the variance of thermal anisotropies in systems in equilibrium on some 
power /3 6 [1, 2) of the volume o 

It has been formally proved 
tern the variance (AMy) 2 of the spatial anisotropies of any extensive magnitude My over 
macroscopically large, but finite, volumes V must grow with their size L as (AMy) 2 ~ L a , 
with a > 2. The derivation of equation (|2T)|) below is a brief proof of this result. Scale 



integration, (AMy) 2 ~ V 13 . 

15^ that in any homogeneous and isotropic statistical sys- 



invariant statistica 
superhomogeneous 



systems like ((H), for which a = 2, have been dubbed for this reason 
16[. We can say that, in this sense, scale invariant systems are the most 
ordered homogeneous and isotropic statistical systems and thermal systems in equilibrium, 
for which a > 3, are more disordered than scale invariant systems. 

Scale invariant anisotropies are known in mathematics, for example, in the distribution 
over the complex plane of zeroes of statistical analytic functions Q|. In physics, scale invari- 
ant anisotropies happen generically in homogeneous and isotropic statistical systems with 
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very long relaxation times to thermal equilibrium. The reader familiar with the arguments 
that justify Eq. Q will appreciate the following illustrative example. 

A. A simple example of scale invariant anisotropies. 

Let consider an homogeneous and isotropic distribution of a large number of individuals 
over a very large (infinite) domain C R 3 and assume that this initial distribution of 
individuals does not contain random density fluctuations with wavelength longer than the 
shortest accessible resolution scale d. 

We now let each individual to walk randomly and assume the same statistical features for 
the individual random walks everywhere in f2. The typical distance D(t) = (\r(t) — rol 2 ) 1 ^ 2 
between the actual position of each individual at time t with respect to its original position 
grows monotonically with the time interval D{t) = ^/j t x l 2 as a function determined by the 
specific properties of the stochastic process. Let us call the distance D(t) a particle horizon: 
at any instant each individual is constrained to be located within a finite region of typical 
size D{t) around its original position. Let us also assume that we have waited enough, such 
that D(t) > d. 

Let us denote by Ny the number of individuals contained within a macroscopic volume 
V C Q at time t. Macroscopic volumes are defined as volumes of typical size L > d. 
Random fluctuations in the extensive statistical magnitude Ny result from the individual 
random walks. We wish to characterize their variance (ANy) 2 as a function of the size L of 
the considered macroscopic volume. 

If the macroscopic volume has a typical size L in the range d ^ L < D(t) (of the order of 
the particle horizon D(t) or shorter) the fluctuations are obviously thermal (AiVy) 2 oc V, as 
we have nothing but an example of brownian motion of independent particles in a macro- 
canonical ensemble. Indeed the particle horizon D(t) is not evident locally, over volumes 
of smaller size. The temperature of these thermal anisotropies is dictated by the statistical 
features of the random walks. If the laws of the individual random walks are invariant under 
spatial translations and rotations in Q the system remains statistically homogeneous and 
isotropic at any time. The temperature of the local anisotropies may change with time, 
but at any instant it is the same over the whole space Q (even though the existence of a 
finite particle horizon implies that the statistical system cannot have relaxed to thermal 
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equilibrium over scales beyond it). 

The existence of a finite particle horizon is manifest in the fluctuations in the number of 
individuals Ny in volumes V of size L ^> D(t). The density anisotropics over volumes much 
larger than the particle horizon are scale invariant (AAV) 2 ~ S, similar to the primordial 
cosmological density anisotropies at decoupling that are the focus of our interest in this 
paper. In order to justify this claim let us first define the cover of width D(t) of the boundary 
surface of a macroscopic finite volume V, and denote it as do(t)V, as the union of all balls 
of radius D(t) around any point x G dV: the cover dr>( t )V of a volume V is the set of all 
points whose distance to its boundary surface dV is smaller than D(t). Individuals that are 
originally located out of this cover cannot contribute at the instant t to fluctuations in the 
number of individuals within the volume V, even though they also are walking randomly. 
Only individuals that were originally located within do(t)V can produce, when they cross the 
border dV, fluctuations in the total number of individuals in V. That is, only the degrees 
of freedom in do(t)V are apparent in these fluctuations. Naturally, then, the anisotropies 
that the random walks of these degrees of freedom induce in the number of particles Ny are 
scale invariant (AAV) 2 ~ D(t) S ~ D(t) L 2 , rather than thermal (AAV) 2 ~ V ~ L 3 . 

The power spectrum of the number density anisotropies generated by the random walks 



of individual particles has been described in the literature 



C, for kZDU)- 1 

V{k-t) = { ' W (17) 

( C[k D(t)]\ for fc<D(t)- 1 , 

where C is a constant determined by the specific features of the random walks. This power 
spectrum is roughly constant over the range of modes with wavelength shorter than the 
particle horizon D(t), but it decreases as ~ k 4 for modes with longer wavelength. The 
suppression of this power spectrum in modes with wavelength longer than the horizon D(t) 
is similar to the suppression of the power spectrum (0) of cosmological density anisotropies 
left at decoupling by standard physics operating during the radiation dominated expansion 
of the universe. Hence, on scales L ^> D(t) condition (|T2j) holds for power spectrum (|17jl. 
with H~q replaced by D(t). In fact, this condition defines the presence of a causal horizon in 
homogeneous and isotropic statistical systems and implies that the system has not reached 
thermal equilibrium over scales beyond it. 

It is clear from this example that the scale invariance of spatial random anisotropies 
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over volumes of size L ^> D(t) is a direct consequence of the fact that individual degrees 
of freedom are constrained within a finite particle horizon D(t). The horizon prevents the 
development of random density fluctuations with longer wavelength. This constrain is ex- 
pressed through the suppression of the power spectrum ()17j) in modes k <C Dit)^ 1 . As a 
result of the strong suppression of the power spectrum in these modes, the largest contribu- 
tion to the integral expression (J3J) for the variance of density anisotropics over volumes of 
size L^> D(t) comes from modes k > _D(£) -1 . Their contribution is scale invariant because 
it corresponds to density fluctuations at the boundary surface of the considered volume of 
integration. 

We wish to remark that the contribution from modes k Z D(t)~ l to macroscopic density 
anisotropies over volumes of size L 3> D(t) is not a mathematical pathology associated to a 
boundary surface that needs to be resolved with infinite precision, as sometimes claimed in 
the literature in cosmology. Scale invariant anisotropies are produced by boundary fluctua- 
tions with typical wavelength d ^ X ^ D(t), much shorter than the size L of the considered 
volume of integration but still much longer than the resolution scale d of the boundary 
surface. 

B. The formal arguments. 

We now wish to formalize these arguments using simple tools and consider an statistical 
density field on R 3 , 

/d^ k ■* 
-—- 5 n {d) e"^, (18) 

to describe, in comoving coordinates, the cosmological density field at the instant of decou- 
pling resolved at a certain finite length scale d. The scale factor of the universe at decou- 
pling a(t eq ) = 1 is chosen as reference. The density field is assumed to be homogeneous and 
isotropic, 

(p(x; d)) = peg (19) 

(or equivalently, (%(d)) = 0), up to tiny statistical fluctuations. We assume that the 
random modes SUd) = S_^(d) are independent statistical complex variables with gaussian 
distribution, 
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(5i(d)% 2 (d)) = (2vr) 3 d) 5 3 (h - k 2 ). (20) 

The function V(k; d) is known as the power spectrum of the statistical fluctuations and is 
always positive V(k; d) > 0. In isotropic statistical systems the power spectrum is a function 
V(k; d) = V(k; d) only of the modulus of the wave vector. 

The finite resolution scale of the density field is introduced in the description through 
a convolution kernel (window function) with typical width d > 0, p(x; d) = J d 3 y W(x — 
y; d) p{y) (see section VIII in |3|). The specific features of the window function W(x — y; d) 
are not of particular interest here, other than it introduces a natural ultraviolet cutoff in the 
the power spectrum V(k; d) = V(k)\W(k; d)\ 2 over scales shorter than the resolution length: 
|W(fc; d)\ 2 <C 1, for k 3> d _1 and |W(&; d)\ 2 ~ 1, for k < d" 1 . In other words, the convolution 
kernel in the definition of p(x; d) simply integrates out all modes with wavelength shorter 
than the resolution scale. 

Macroscopic volumes are defined as volumes whose typical size is much larger than the 
resolution scale d. According to the central limit theorem the statistical extensive magni- 
tudes 

M v (d) = f d 3 x p(x;d) (21) 

have gaussian distribution. Hence, they are completely specified by their first two momenta: 
the average value (My) = p eq V and the variance of its random fluctuations (AMy) 2 , which 
is related to their power spectrum by expression (JSJ). 

Expression (JHJ) emphasizes that the variance of random density anisotropies over macro- 
scopic volumes of any size L receives positive contributions from random fluctuations in all 
modes <%(c£). We want to evaluate separately the contribution from modes k 3> Lr x and 
compare it to the contribution from modes k < L^ 1 , 
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(22) 



The first contribution, from modes k < 1/L, coincides with approximation ©: 

rk<L~i d 3 k 
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\-T(k ~ 1/L) K 2 ~ ~ 1/L) V. (23) 
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The second contribution, from modes k ^> L 1 , is neglected in approximation This 
contribution is the focus of our interest in this paper because it is scale invariant, 



k<d" 



d 3 k 
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d x e 



-ik-x 
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fc»L- 1 fc 4 



5. 



(24) 



As above, S stands for the area of the surface that bounds the volume of integration V. 
Equation (|2*4*j) is very simple to derive in the case when the volume of integration V is a 
sphere of radius L, 



d 3 x e 



-ik-x 



2 S f 1 

= Air — I - — sin(kL) — cos(kL) 
fc 4 \kL 



(25) 



and the power spectrum V(k) in the range k ^> L~ x is smooth over scales dk < L^ 1 . 
The result can be immediately extended to any other connected volume, assuming that the 
surface that bounds it is smooth and its area does not depend on the resolution scale. 

In summary, the variance of spatial random anisotropies of extensive magnitudes over 
macroscopic volumes V of any size 



(AM V 



p 2 V(k ~ 1/L) V + 



Plq 
7i 



k 2 



S. 



(26) 



receives always a bulk contribution from fluctuations with wavelength comparable to the size 
of the volume of integration and an additional scale invariant contribution from fluctuations 
at its boundary surface with much shorter wavelength. We notice that the global amplitudes 
of each of these contributions depend, respectively, on the power spectrum in modes k < 1/L 
and k 3> 1/L. Obviously, if the first contribution is negligible compared to the second, then 
the anisotropies are scale invariant. 

Let us now consider, in particular, the density anisotropies produced at decoupling by 
standard physics operating during the radiation dominated expansion of the universe. Their 
power spectrum @ is strongly suppressed in modes k H eq with comoving wavelength 
much larger than the causal horizon. Over cosmologically large comoving volumes, of size 
L ^> H^ 1 , condition (jT2j) holds. An explicit evaluation of the two contributions in (J2lj|) gives 



p% V(k ~ 1/L) V < P % V(k ~ H eq ) 
for the bulk contribution and 



V ~ p% V(k ~ H eq ) H~ 4 L- 1 



eq) - t - t eij 



(27) 
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S * p% V(k ~ H eq )H; q l L\ (28) 

for the scale invariant contribution. We notice that ()27|) corresponds to the estimation (|1U|) 
of the variance of primordial cosmological density anisotropies at decoupling in standard 
FRW cosmology, which led to the formulation of the origin of structures problem. This 
estimation neglects the second contribution (J28|) . Nevertheless, a direct comparison clearly 
shows that over volumes of comoving size L 3> H~ q the second, scale invariant contribution 
(|2*5J) . from modes k Z H eq , prevails: 

(AM V ) 2 ~ p% V{k ~ H eq )H^L\ (29) 

while the first, bulk contribution (|27jl . from modes k < 1/L, is negligible. The flawed 
estimation (|1U|) keeps the negligible contribution (J27)) and neglects the dominant one (J28|) . 

The power spectrum V(k ~ -f^eg) i n modes with comoving wavelength comparable to the 
comoving causal horizon is not suppressed with respect to the power in modes with shorter 
wavelength. Therefore, over comoving volumes V eq of the size of the horizon approximation 
© is valid, (AM Veq ) 2 ~ p 2 eq V(k ~ H eq )V eq ~ p^P(A; ~ F e(? ) F" 3 . Equation ((131) is obtained 
straightforward from (|29j) . 

Finally, we wish to remark that in the framework we have laid in this section we did 
not even need to worry about the technical difficulties that appear when trying to manage 
the coupled dynamics of density fluctuations with comoving wavelength longer than the 
causal horizon and space-time metric perturbations because the power spectrum in 
these modes is strongly suppressed and their contribution to cosmological anisotropies is 
sub-leading. 

C. Scale invariant gravitational potential. 

We have claimed in the Introduction to this paper that scale invariant primordial den- 
sity anisotropies over cosmologically large comoving volumes (I1I2J) produce a statistically 
constant gravitational potential over cosmologically large comoving distances. In this 
subsection we prove that this result is a consequence of the strong suppression (jT2*j) of the 



Peq 
71 



dk 



V(k) 



s 



Ti 



dk 



P(k) 
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power spectrum of primordial density anisotropies in modes with cosmologically large co- 
moving wavelength. 

The statistical gravitational potential at the time of decoupling 

e~ lk -\ (30) 



(2tt) 3 Yk 

is related to the statistical density field at the time (|TH|) by the equation 

V 2 0O?) =4irG(p(x)- Peq ). (31) 
In momentum space the same equation can be written 

cj> % = -4nG Peq -L %. (32) 



The gravitational field (J3*U|) is real because 

P-% = fa 

and it is zero everywhere on average (<f){x)) = 0, because 



(33) 



(fa) = 0- (34) 
Thus, the average potential difference between any two points is also zero, 

(0(f) - 0(y)> = 0. (35) 

The two points correlation function of the statistical gravitational potential is related to the 
power spectrum of density fluctuations by the expression, 



(2vr)3 7 (2tt)3 1^121^ 



(27r) 3 7 (2tt) j 

=(^) a /S¥ cH ^ M (36) 

and the variance of random fluctuations in the potential difference between any two points 
is given by, 
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([0(f) - <p{y)f) = 2(0(f) 2 ) - (m m) - {m = 

We are interested in the asymptotic behaviour of this variance when \x — y\ ^> H~ q) the 
comoving distance between the two considered points is cosmologically large. To evaluate the 
integral expression (J37)) we separate the contribution from modes k <C H eq , whose comoving 
wavelength is much larger than the horizon, from the contribution of modes k > H eq , whose 
comoving wavelength is within the horizon, 



f)-0(y)] 2 )^(4Gp, 



eqj 



Jo k \ k\x — v\ / Jk>H K „ K 



(38) 



Jo k 2 \ k\x — y\ J Jk>H eq 

The latter contribution is roughly constant {AGp eq ) 2 p(A ^ eg) because Sin k \ k J^ l] < 1, for 
k > H eq and \x — y\ ^> H~ q . We want now to explore the contribution from comoving modes 
with cosmologically large wavelength when the power spectrum of density fluctuations is 
largely suppressed in this range (fT2J). 



rk« Heq dk THftf_ Sin [k\x-y\] \ < V(k ~ H eq ) r*<*.« ^ ^ f _ Sin [k\x - y\] \ 
Jo k 2 \ k\x — y\ ) ~ Jo \ k\x — y\ J 

T(k ~ H eq ) 



< 



H eq 



We find again that the contribution from random density fluctuations with cosmologically 
large comoving wavelength is suppressed relative to the contribution from density fluctua- 
tions with comoving wavelength wihin the horizon. In consequence, (<f>(x)<f>(y)) "C ([4>{x)} 2 ) 
when \x — y\ ^> H~ q and 



*o - mf) * (W)?) + {[m?) * ^G Peq y p( ^ geg) (39) 

the random gravitational potential is statistically constant over cosmologically large comov- 
ing distances. 
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D. Scale invariant anisotropies vs. thermal anisotropies. 

The result ()13p that we have proved in subsection II. B above is at odds with the intuition 
inferred from the common experience with thermodynamic systems in equilibrium, in which 
the variance of random density anisotropies of extensive magnitudes grows with the size of 
the considered region of integration as some power of its volume (AMy)t hermal ~ V 13 , with 
1 < (3 < 2. In this section we want to compare the bulk dependence of the variance of thermal 
anisotropies with the surface dependence of the variance of scale invariant anisotropies. 

Consider a thermodynamic system in equilibrium far from a phase transition. Thermal 
fluctuations in the local densities are described by Poisson two points correlation functions 
F(x — y) = \ [(p(x)p(y)) — (p(x)) (p(y))] = ( 2 d 3 5 3 (x — y), or the more realistic 

F(r) = C ^ e-^ rf ) 2 , (40) 

where d ~ 1/T is the finite resolution scale and ( is a dimensionless factor that measures 
the global amplitude of the fluctuations. The corresponding power spectrum, 

V(k) = J d 3 rF(r) e +i ^ ~ C 2 d 3 (41) 

is roughly constant over the whole range < k < d~ l of accessible scales. Hence, over 
spheres of radius L the variance of spatial thermal anisotropies is 



r L / d I 1 1 \ 

(AM,) 2 = 8C 2 pi d 3 L 3 J Q d{kL) \—-sin{kL) - ^cos(kL)j ~ C 2 Po d 3 V. (42) 

It can be readily seen that, independently of the size of the region L, the largest contribution 
to the last integral comes from modes kL ~ 1 (the integrand vanishes for fcl C 1 and for 
kL ^> 1) and the variance is linearly proportional to the volume of integration, in good 
agreement with the estimation (jUJ). 

The linear dependence of the variance ()42|) on the volume of the considered region of 
integration results from the short range of the Poisson two points correlation function ()40j). 
Poisson correlation is characteristic of thermodynamics systems in equilibrium far from a 
phase transition. For thermodynamic systems in equilibrium but close to a phase transition 
the power spectrum of thermal Actuations is peaked at modes with very long wavelength 



20 



V(k) ~ C 2 d 3 7 k 7 , with 7 e (0,3), and the associated two-points correlation function is 
long-range. When this is the case 

r L / d \ ( \ 1 \ 2 

(AM,) 2 = 8« ^ / d(kL) ^ [ m ^n(kL) - -co a{ kL)} 

~ C pi d 3 -' V 1 *-"' 3 . (43) 

The largest contribution to the variance of spatial density anisotropics still comes from 
modes kL ~ 1 whose wavelength is comparable to the size of the volume of integration, but 
now the variance grows with a higher power of the volume of integration. Estimation © 
still gets it right. 

We have brought these examples of thermal anisotropics into consideration because we 
want to notice that their power spectra are not suppressed over any range of short momenta, 
limk-,oP{k) ^ 0. In the statistical systems described by these power spectra there is no 
causal horizon. These systems are theoretically allowed to excite with finite probability 
random density fluctuations with any wavelength longer than the resolution scale d. In 
consequence, in these statistical systems the bulk contribution (|2*7j) to the variance (J2T?|) of 
spatial random anisotropies in macroscopic volumes V of any size L is always dominant over 
the scale invariant contribution (|28p. 

The power spectrum of an statistical system that is constrained by a finite causal horizon 
D looks like 

V(k) = p(k D) V (k ~ D^ 1 ), (44) 

where the dimensionless function p(k D) ~ o(k D) decreases to zero faster than linearly for 
k D^ 1 and p(k D) ~ 1 for k ~ £) -1 . The pattern of anisotropies associated to these 
power spectra are typically of the form fi2|) or P3|) over volumes of size L < D shorter than 
the horizon, but it turns scale invariant (fT3*|) over volumes of size L > D as a result of the 
suppression of the power spectrum ()44|) in modes k <C D~ x 

In the model discussed in subsection II. A, for example, the power spectrum ()17|) is sup- 
pressed in modes k ^ D^ 1 by a factor p(k D) ~ (k D) 4 . A similar suppression factor 
appears in the power spectrum of cosmological density anisotropies produced at decoupling 
by standard physics operating during the radiation dominated expansion. A suppression 
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factor p(x) = jt^t, n > 1, appears in the power spectrum that describes the statistical 
distribution of topological defects in a sample cooled down homogeneously, but very fast, 
below its critical temperature (Kibble mechanism) Q|. The super fast cooling rate of the 
sample, rather than the cosmological expansion, introduces an "horizon" in this example, 
if we define the horizon D of an homogeneous and isotropic statistical system as the scale 
D^ 1 beyond which its power spectrum is suppressed. Spatial random anisotropies in the 
distribution of defects over regions of the sample of size much larger than the "cooling" 
horizon are theoretically known to be scale invariant. 

The suppression of the power spectrum (J4*lj) over scales k <C D^ 1 introduces the global 
constrain 



This condition f Q d 3 r F(r) ~ 0, which we have associated directly to the scale invariance 
of random anisotropies over very large volumes, is known in the literature as a demand of 
incompressibility of the statistical density fluctuations. Incompressibility is understood in 
our context as a constrain (AMq) 2 = of total conservation of the extensive magnitude M 
over the whole volume Q, see section IV in Q|. This condition can be evidently satisfied 
when the sub-volumes V of any size are assumed to be embedded in a larger isolated micro- 
canonical ensemble. In the framework of a quantum field theory this condition can be easily 
obtained if the state of the system is assumed to be an eigenstate of Mq. In quantum Hall 
effect in condensed matter physics, for example, the incompressibilty of the electron fluid is 
known to be associated to the low-energy border exitations of the system 



In this subsection we examine the amplitude of the scale invariant cosmological density 
anisotropies left at decoupling by standard physics (fl~3|) and compare it with the amplitude 
of the scale invariant anisotropies (0) imprinted in the CMBR. 

We describe the radiation dominated density field at (physical) cosmological time t < t eq 
as a statistically homogeneous and isotropic random field, 





(45) 



E. Amplitude of primordial cosmological density anisotropies. 
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p(x ; t)=pt + p t 



d 3 k 

(2^3 V e 



■ik-x 



(46) 



with 



<M = ' (M*> *)> = *), ( 4? ) 

and 

(<^ t .£ 2 > = (2tt) 3 P(N;t) - fc 2 ). (48) 

The Fourier modes are assumed to be statistically independent and gaussian. We keep using 
comoving spatial coordinates defined with the scale factor of the universe at decoupling 
a(teq) = 1 as reference. 

The statistical variable p(x; t) — p(y; t) describes the random fluctuations in the density at 
some comoving site x with respect to a different comoving site y, at the given cosmological 
time t. Statistical homogeneity (|47|) demands that this statistical variable has zero average 
value, 

(p(x;t)-p(y;t)) = 0. (49) 
We are interested in the variance of its statistical fluctuations, 



{[p(x;t) - p(y;t)} 2 ) 



([p(x;t)} 2 ) - (p{x;t)p{y;t)) 



2pt 



([p(x;t)} 2 ) ~ Pt 
d 3 k 



V(k;t) (l 



— e 



(p{x; t)p{y;t)) - p\ 

+ik-(x—y) 



71^ 



q- dkk 2 V(k;t) 1 



Sin(\k\\x — y\) 



\k\\x 



y\ 



(50) 



During the radiation dominated expansion of the universe the cosmic fluid supports the 
propagation of non-linear density waves, which travel with sound velocity c s = 1/ (in 
natural units) and operate to erase density inhomogeneities. Therefore, we must expect 
maximal statistical correlation between causally connected comoving sites: 
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<[p(£;t)-p(y;t)] 2 >«2 ([p{x ; t)} 2 ) - p 2 



for \x — y\ < H t . 



(51) 



On the other side, causality constrains forbide the existence of statistical correlation between 
causally disconnected sites: 



(\p(x; t) - p(y; t)f) ~ 2 [<[p(f; t)} 2 ) - p 2 t \ , for \x - y\ » H' 1 . (52) 

Actually, causality does constrain the two points correlation function to satisfy the more 
stringent condition ()45j) . The comoving causal horizon Hf = (t/teq) 1 ^ 2 grows mono- 
tonically during the radiation dominated expansion of the universe. 

These features correspond to a power spectrum V(k; t) peaked at k ~ H t , 

V(k~H t ;t)=P t , (53) 
and suppressed both at k ^> H t and at k -C H t , 

, V t Hf/k 3 , for k > H t , 
V{k;t)<{ 11 ' 1 ' (54) 

V t k 4 /Hf, for k < H t , 

We plug this power spectrum in expression (J5J) to obtain the variance of the density 
anisotropies left in the cosmic plasma at cosmological time t: 

I ) The suppression of the power in modes k ^> H t implies that over volumes V of 
comoving size L < Hf smaller than the causal horizon the anisotropies are extensive, 

(AM v ) 2 (t) ~ p\ (H?V t ) V 2 ~ (AM Vt ) 2 ^ 2 . (55) 
In particular, over comoving volumes Vt of size L ~ H^ 1 comparable to the causal horizon 



(AM Vt ) 2 (t) ~ (p^ t ) V 2 ~ p^r 3 ^- (56) 

/J) The suppression of the power in modes k H t implies that over volumes V of 
comoving size L ^> if^r 1 much larger than the causal horizon the anisotropies are scale 
invariant 
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(AM v ) 2 (t) ~ p 2 /^ L 2 ~ (AM Vt f-^. (57) 

The amplitude of the cosmological density anisotropies (jB3j) , (J57)l at decoupling t = t eq is 
commensurate with the amplitude of the primordial anisotropies imprinted in the CMBR if, 
and only if, condition (fTljl is fulfilled. From (JHHJ) we obtain the condition H% Vt eq ~ 10 -ltl 
In order to give it a clear physical meaning we now define 



2_ 1 



(ip(x ; t)] 2 ) - p 2 ] = / ^3 «) ^ A? Vu (58) 

and notice that it is the normalized variance of the density fluctuations between comoving 
sites x and y that are causally disconnected at time t, 



K 



(\p(x; t) - p(y; t)] 2 ), for \x - y\ » H; 1 . (59) 



1 1 

We have from (|56[) and 



(AM Vt ) 2 (t) ~ n 2 t p 2 t V t 2 . (60) 
Hence, condition (|14jl simply requires 

< « 10- 10 . (61) 

As there does not seem to exist any obvious theoretical or observational constraint that 
rules out (J6*T]) in the framework of standard FRW cosmology, we should conclude from 
this analysis that primordial cosmological structures $1} could have been seeded during the 
radiation dominated expansion of the universe. 

In a forthcoming paper we will explore the mechanisms involved in actually fixing n t at 
the time t = t eq of decoupling. We will show that k 2 is a measure of the entropy of the 
random density fluctuations in the cosmic plasma at time t. This observation is important if 
the entropy of the density fluctuations is preserved during the adiabatic radiation dominated 
expansion of the universe. In such constant since the much earlier time i* when 

the universe entered the stage of adiabatic expansion, 



k 2 — constant, for t G (t*,i e g)- (62) 
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This would imply that the amplitude of the anisotropies (|55|). (|3B)l. (jBTf) at decoupling 
t = t e(? was indeed fixed by entropy-producing mechanisms at that earlier time t*. Notice 
that condition 1)62)1 is satisfied if the relative amplitude of the density fluctuations between 
any two causally disconnected sites (|HHJ) is preserved during the adiabatic expansion. 



We have shown that the prevalent analysis of the origin of primordial structures in stan- 
dard FRW cosmology is flawed by an unjustified approximation, which has led researchers 
to erroneous conclusions. Our revised analysis of the issue shows that standard physics 
operating in the cosmic plasma during the radiation dominated expansion of the universe 
can produce scale invariant cosmological density anisotropies at decoupling with amplitude 
comparable to that of the primordial anisotropies imprinted in the CMBR. Indeed, the char- 
acteristic scale invariance of the primordial cosmological density anisotropies at decoupling 
seems to be a signature of the causal structure of the standard FRW universe, which prevents 
the development in the radiation dominated cosmic plasma of large density fluctuations with 
cosmologically large comoving wavelength. This might explain the reported absence of sta- 
tistical correlations at large angular scales in the temperature anisotropies measured by the 
WMAP hfl- 
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